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Abstract. The relation between saddle points of the potential of a classical 
many-particle system and the analyticity properties of its Boltzmann entropy is 
studied. For finite systems, each saddle point is found to cause a nonanalyticity in 
the Boltzmann entropy, and the functional form of this nonanalytic term is derived 
for the generic case of potentials having the Morse property. With increasing 
system size the order of the nonanalytic term grows unboundedly, leading to 
an increasing differentiability of the entropy. Nonetheless, a distribution of an 
unboundedly growing number of saddle points may cause a phase transition in 
the thermodynamic limit. Analyzing the contribution of the saddle points to the 
density of states in the thermodynamic limit, conditions on the distribution of 
saddle points and their curvatures are derived which are necessary for a phase 
transition to occur. With these results, the puzzling absence of topological 
signatures in the spherical model is elucidated. As further applications, the phase 
transitions of the mean-field XY model and the mean-field fc-trigonometric model 
are shown to be induced by saddle points of vanishing curvature. 

Keywords: Classical phase transitions (theory), energy landscapes (theory), solvable 
lattice models 

1. Introduction 

Phase transitions are abrupt changes of the macroscopic properties of many-particle 
systems under variation of a control parameter. Typical examples are the sudden 
disappearance of the electric resistance when cooling a superconducting material 
below its transition temperature, or the evaporation of a liquid at temperatures 
above its boiling point. An approach commonly used for the theoretical study of 
phase transitions is the investigation of the analyticity properties of thermodynamic 
functions like the canonical free energy of enthalpy. It is long known that non- 
analytic behaviour in a canonical thermodynamic function can occur only in the 
thermodynamic limit in which the number of degrees of freedom N of the system 
goes to infinity [Ij. 

Many researchers took it for granted that the same were true also for micro- 
canonical thermodynamic functions. Recently, however, it was observed that the 
microcanonical entropy, or Boltzmann entropy, of a finite system is not necessarily 
real-analytic, i.e., not necessarily infinitely many times differentiable. Nonanalytic 
entropy functions of finite systems have been reported for certain classical models 
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[21 [31 m [S] as well as for quantum systems |6J, where the latter result relies on a 
suitable but rather unconventional definition of the density of states. 

In light of their conceptual similarity to the definition of phase transitions, it is 
tempting to regard finite-system nonanalyticities of the entropy as phase transition-like 
phenomena. This point of view is advocated in [3] , and the authors of that reference 
argue that such nonanalyticities should also be measurable experimentally, at least in 
very small systems. However, such an interpretation is complicated by the fact that, 
as discussed in [5] , for typical models the number of nonanalytic points of the entropy 
increases unboundedly with the number of degrees of freedom N. 

Due to their typically large number, one might assume that nonanalyticities of 
the finite-system entropy were unrelated to the occurrence of a phase transition in the 
thermodynamic limit. From a theorem by Franzosi and Pettini [7], however, it can be 
deduced that — at least for a certain class of short-range systems — a relation between 
finite-system and infinite-system nonanalyticities exists (and we will come back to that 
theorem later). The purpose of the present article is to further clarify and quantify 
the relation between saddle points of the potential energy and nonanalyticities of 
the entropy for the cases of finite as well as infinite systems. For infinite systems, 
at least in the common situation of ensemble equivalence, a nonanalyticity of the 
entropy will correspond to a nonanalyticity in the canonical free energy, and hence to 
a phase transition according to the standard definition. Note that in the literature 
other approaches relating saddle points of the potential to phase transitions have been 
proposed. Examples include the work of Rutkevich [8], Wales [9], and Pettini [TO] . 
and we will elaborate on the relation to the latter approach in Sec. [H 

Our approach is the following: After fixing notations in Sec. [21 we discuss the 
analyticity properties of the Boltzmann entropy for a generic class of classical many- 
particle systems in Sec. [3] For the case of finite systems, our main result is that 
every saddle point (including maxima or minima) of the high-dimensional potential 
energy landscape corresponds to a nonanalyticity of the Boltzmann entropy, and the 
functional form of this nonanalytic term is derived explicitly. (A weaker version of 
this result has been announced without proof in a recent Letter [H]). This result 
gives a model-independent and quantitative account of nonanalyticities in the finite- 
system entropy as observed for the special cases in [2l [3l [5]. Then, the density of 
states is split into two terms: The first is a sum of the nonanalytic contributions 
stemming from the saddle points, whereas the second contains the analytic rest. 
Performing the thermodynamic limit of the thus obtained expressions, we obtain in 
Sec. [H a condition on the distribution of saddle points and their curvatures which 
necessarily has to be fulfilled in order cause a phase transition. (This result has been 
outlined in the recent Letter [12] ■) We apply our findings to the spherical model 
of a ferromagnet in Sec. [3 explaining the puzzling behaviour that, for this model, 
the topological signatures observed do not coincide with the phase transition energy 
[131 [13]. As further examples, the mean-field XY model (SecHj) and the mean-field k- 
trigonometric model (Sec. [7]) are discussed, illustrating nicely the relation of curvature 
properties of saddle points of the potential and the occurrence of phase transitions in 
the thermodynamic limit. We summarize our findings in Sec. [H 
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2. Preliminaries 

We consider classical systems of N degrees of freedom, characterized by Hamiltonian 
functions of standard form, 

1 ^ 

H{p,q) = -J2p^ + V{q), (1) 
1=1 

where p — (pi, . . . ,pn) is the vector of momenta and q — (gi, . . . , (/at) the vector 
of position coordinates. The restriction to a quadratic form in the momenta is not 
essential, but simplifies the presentation. The potential V is an analytic mapping from 
the configuration space Fat C onto the reals. Whenever Fjv is noncompact, we 
assume V to be confining, i.e., 

lim V{Xq) = oo VO 7^ g G Tat. (2) 

A — >co 

In general V will have a number of critical points (or saddle points) qc, defined 
as points from Fat with vanishing differential, dV{qc) = 0. Note that throughout this 
article the expression "saddle points" is used synonymously to "critical points" and 
includes also minima and maxima. We will assume in the following that all critical 
points of V are non-degenerate, i. e., that the determinant of the Hessian f)v fulfills 

det[^v{qc)]^^ VgcOfF. (3) 

In this case, V is called a Morse Junction. Conceptually, this is an insignificant 
restriction, since Morse functions on some manifold M form an open dense subset in 
the space of smooth functions on M jTSl H^. Therefore, if the potential V of the 
Hamiltonian system is not a Morse function, we can deform it into a Morse function 
V by adding an arbitrarily small perturbation, e.g. 

N 

V{q)^V{q) + Y.^,q, (4) 

i=l 

with small /ij £ IR (i = 1, . . . , iV). For practical purposes, adding a perturbation — and 
thereby destroying a symmetry present in V — may, however, render the computation 
of critical points and indices more complicated or even impossible. 

In case V is not a Morse function due to the presence of a symmetry, an alternative 
strategy is possible: zero-eigenvalues of the Hessian caused by, say, translation 
invariance can be removed by fixing one or several position coordinates, thereby 
destroying the symmetry. Thermodynamic functions of such a system with a reduced 
number of degrees of freedom will differ from the original ones only by a physically 
irrelevant additive constant. 

3. Finite system nonanalyticities 

When computing thermodynamic properties of a Hamiltonian system, a quadratic 
form in the momenta as in the Hamiltonian ([T]) merely leads to a shift in the free energy. 
Hence, for the discussion of phase transitions or nonanalyticities of thermodynamic 
function we may disregard the kinetic term and concentrate on the configurational 
part (see [T7] for a discussion of the pitfalls of this reasoning). 
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Our aim is to investigate the analyticity properties of the configurational density 
of states as a function of the potential energy per degree of freedom v = V/N, 

nN{v) = / dqS[V{q)-Nv], (5) 

or, equivalently, of the entropy 

SN{v) = j^\n[nN{v)] (6) 

of a classical system of N degrees of freedom, characterized by a Morse function V. 
The crucial observation for our analysis is that nonanalyticities of rijv are closely 
related to the critical levels of V, i.e., the values Vc of the potential per degree of 
freedom V/N for which a critical point qc exists such that V{qc)/N = Vc- If [vi, W2] is 
an interval free of critical levels, one can show that and are smooth functions 
in this interval (Lemma 2 of TS]). 

The Morse property ([3|) guarantees that all critical points of V are isolated. We 
thus may reduce the discussion to the effect of a single critical point qc of V on the 
analyticity properties of the density of states f^Ar. Without loss of generality we can 
choose V{qc) — 0. Furthermore, making use of the Morse lemma [151 IHl) a local 
coordinate system x — {xi, X2, ■ ■ ■ , xn) can be chosen such that 

k N 

V[q{x)] = -Y,^!+ E (7) 

1=1 i=fc+l 

in some open neighbourhood of Qc = q{0)- We denote the Jacobian of the 
transformation into this coordinate system by ^(x) and its determinant by J{x). We 
will use the expansion of J{x) at a; = for the derivation of the singular behaviour of 
flN^v) at u = 0, 

using a multi-index notation The zeroth order / = is particularly 

symmetric and leads to the leading nonanalytic behaviour of VIn- In this section we 
will restrict the discussion to / = 0. The analysis of higher order corrections demands 
different methods and is shifted to the Appendix. The coordinate transformation at 
gc = is linked to the second derivatives of the potential. Thus J(0) can be written 
as 

J(0)=ao = |det[i3y(0)/2]r^/^ (9) 

The determinant det9jv can be viewed as a measure of some "curvature" at a critical 
point, and this allows us to interpret J in our later results as quantifying the flatness 
of a critical point. 

The index k in Eq. ^ is called the Morse index of the critical point. It equals the 
number of negative eigenvalues of the Hessian ^v'(O) of V at the critical point qc ~ 
and it determines whether a critical point is a minimum (fc = 0), a maximum (fc — N), 
or a proper saddle (fcs{l,...,iV — 1}). In this section we will mainly deal with the 
latter case. The simpler calculation for the extrema is entailed in the more general 
set-up of the Appendix. Following the notation in ^20], we introduce new variables 
X > and y > by defining 

k N 

x^ = Y.^l Y'^Y. (10) 

i=l i=k+l 
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V = -X^ + Y^. (11) 

We restrict the calculation of il^ to a ball of radius r centered at x = 0, 

X2 + y2<^r2, (12) 

and we choose r small enough such that the ball fits into the open neighbourhood 
of Morse's lemma||| The restriction to the ball will result in spurious singularities at 
Nv = ±r^ which are the limiting values of V on the ball. However, we know that the 
relevant nonanalyticity has to be located at v — and thus may assume N\v\ < in 
the following. 

Introducing polar coordinates for [xi, ...,Xk) and (xfe+i, xjv), we find for the 
density of states in the vicinity of qc to leading order the expression 

n%^\v) = J(0) / d^x 5 {-X^ + Y^~V)Q (r^ - - Y^) (13) 



= J{0)Ck-iCN-k-i / dX / AYX'^-^Y^-''-^ 
Ja Jq 

X ,5 {-X'^ + Y'^ ~V)e (r^ -X'^ - Y^) , (14) 
where C„ is the volume of the n-sphere, 
27r(»+i)/2 
^" ^ r[(n + l)/2]' 

and 

f if x < 0, 

denotes the Heaviside step function. We evaluate the X-integral making use of 

6{-X^ + Y^ -V)= S[2{Y^ - vy^^{X-[Y^ - V]^/^)]eiY^ - V), (17) 

obtaining 
with 



1 l.y/(r^+V)/2 

i± = - dy(r2_i/)(^-2)/2^^,_fc_i^ ^^g^ 

2 Jvv0(v) 

In order to avoid cuts while comparing the two integrals /„ and /+, we consider V as 
a complex variable and shift it slightly, say, above the real axis, 

ImV" = e > 0. (20) 

We use the variable transformation 



Y = ^Vy (21) 

where we interpret the square root as a single valued function on the complex plane 
cut from — oo to (and \/T = +1). With this definition of the square root we resolve 

I The precise shape of the region is not essential for the calculation. 
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the ambiguities in raising complex numbers to half integer powers. The integrals I± 
are transformed into 

/_ = - y Av^V{y-\) , (22a) 

/+ = J dyVV{y - 1) VV^ . (22&) 

Since ImV^ > 0, the upper limit of the integral lies in the lower half plane The 
integrand is holomorphic in , allowing us to shift the contour freely in the lower half 
plane. The contour is supposed to only touch the real line in the lower limit ?/ = 1 or 
y = 0, respectively. By our convention for the square root we have \/y — 1 = — i^l — y 
in If, in the case of /+, we let the contour take a detour via — ie, we are able 
to express both integrals in terms of incomplete beta functions B(fci, fc2, z), 

-N-2 



V 



4[k- 



2 



' 2 ' 2' 2F j " V~2 '2' 



(23 a) 



VV fN-k k Y'^ + V 



^-=^i^n^'2'^j- ^'''^ 

It is also possible to express the integrals in terms of beta functions with arguments 
in the interval (0, 1). To achieve this for /„, we have to apply a Mobius transform 
y = z/(z^ 1) which permutes the points (0, 1, oo) to (0, oo, 1). For the transformation 
of /+ we have to permute (0, 1, oo) to (oo, 0, 1), facilitated by y = 1/(1 — z). Performing 
these substitutions we obtain 

y(7V-2)/2 fk N r^-V\ , , 

1+ = B -,1 , (24a) 

We obtain for the difference of /+ and /_ 

+ 41*^ r(iV/2) ' ^ ' 

where the complete beta function has been expressed in terms of the F-function. 

What are the implications of this result for the nonanalyticity of I{V) — 
&{V)I+ + 8(— y)/_? We know that I{V) is a real function which may have a square 
root singularity (for odd N) or a logarithmic singularity (for even N) at F = [this 
follows for example from the properties of the beta functions in Eqs. (|24a[) and (j24&P ]. 

In the case of even k, the difference /+ — /_ is real for V > and we simply 
conclude that I{V) is analytic up to a term (-l)'=/2y(^-2)/2e(y)r(fc/2)F((A^ - 
fc)/2)/[4F(7V/2)]. 

In the case of odd k and odd N, we find that — /_ is real for < 0: We have 
to keep in mind that for ImV > we have \^—V = —iVV to see that liV) is analytic 
up to a term {-l)^^-''^/^-V)^^-^'>/^e{-V)T{k/2)T{{N - fc)/2)/[4F(A^/2)]. 

Finally, in the case k odd and N even, — /_ is imaginary for positive and 
negative V. This is a consequence of a logarithmic singularity which has to be present 
in both, /+ and /_ (otherwise Eq. ([25|) would contain an In(V^)-term). In fact, for 
ImV > a term ln|T^| in I{V) results in a term ln|y| + iarg(t/) in /+ and a term 
\n\V\ — i[TT — aig{V)] in /_ which gives in in 1+ ~ Comparison with Eq. (p5)) 
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yields that I{V) is analytic up to a term {-l)'-''+'^'^/'^V'^'^~^^^'^{\n\V\)r{k/2)r{{N ~ 
fc)/2)/[47rr(A^/2)]. We confirm that the nonanalyticity is independent of the radius r 
of the domain under consideration. 

The leading order nonanalyticity depends on the remainder of the index k mod 4 
only. This remarkable fact prevails to all orders as we will show in the Appendix. 
Moreover, the type of the nonanalyticity depends only on whether N and k are even 
or odd. 

Including the prefactors in Eq. (fT4| and going back to v — V/N , we obtain the 
leading order behaviour of the density of states as summarized in the following 
proposition. 

Proposition 1. Let V : G Hi be a Morse function with a single critical point of 
index k in an open region G. Without loss of generality, we assume V{qc) = 0. Then 
the density of states riN{v) is nonanalytic atv = 0. The leading order nonanalyticity 
is given by a function depending on the number of degrees of freedom N, the index of 
the critical point /cmod4, and the potential energy per particle v. We have 

"^'^"^ ^ ivr(iv/2V|dL[iov-(gc)/2]|''""''=-°'^^^"^ ^''^ 

with the universal function 

r (_l)fe/2^(A'-2)/20(^) for fc even, 

[ (_l)(W-fc)/2(_„)(A^-2)/20(_^) ^ 

In fact, it is possible to give a complete description of the nonanalytic part of the 
density of states in terms of the Taylor coefficients of the Jacobian determinant ^ 
at the critical value. We find that only those indices / — {ii, . . . , zat} contribute to 

where all ik (fc = 1, . . . , TV) are even. The general situation is summarized in the 
following theorem. 

Theorem 1. Let V : G M, be a Morse function with a single critical point Qc of 
index k in an open region G and V{qc) = 0. The density of states can be decomposed 
into an analytic part Qff and a nonanalytic part fl'}^, 

The nonanalytic part of the density of states is the product of an analytic function 
^Nkj '^'^^ universal nonanalytic function /ijvfcmod4 ff'^s*^ ^Q- j^Tp - The 
analytic factor depends on the number of degrees of freedom N , the index of the critical 
point k, the Jacobian determinant J of the transformation from the original coordinate 
system to Morse coordinates, and the potential energy per particle v. The nonanalytic 
factor depends on N , k mod 4, o,nd v, 

^N{^)^h-^.kA^)h'}^,kraoUv)- (29) 

The analytic factor /i^ j{v) can be expressed in terms of the expansion coefficients 
ai of the Jacobian determinant at the critical value, Eqs. 0^ and With Li ~ 

ii + . . . + ik, L2 = ik+i + . • . + In, and \L\ = Ii + J2 we have 

^W/2^W/2+|/|-lJ|A^ (^^, 
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j 1 if or fc even, 
^ ~ j 2 if TV and fc odd. ^^^-^ 

Proof. Eq. ([50]) will be derived in the Appendix. The leading order / = yields the 
prefactor in Eq. (p6)) . thus proving Proposition [T] in the cases fc = and k — N. It 
remains to show that j, j{v) is analytic at w = 0. To this end, it is sufficient to show 
that the sum in Eq. pO|) converges in some neighbourhood of w = 0. We may assume 
that the Jacobian determinant is analytic in a neighbourhood U of x = 0, yielding 
|a/| < where R is the maximum of the Xi in U. Using Sterling's formula we find 

that, in the limit |/| — > oo, the summands in Eq. I|30p are bounded by 

yielding li^/N for the radius of convergence at u = 0. □ 

The theorem gives a complete account of the nonanalyticities in Ojv as the 
consequence of a single critical point of a Morse function V. In the presence of several 
isolated critical points, their nonanalytic contributions simply have to be added up. It 
is also possible to state the result in terms of the potential V in the original variables, 
not making use of Morse coordinates. The result, however, is more complicated and 
for our purposes the above version will be sufficient. 

One can verify that in any of the three cases in (|27|) . f^Ar is [{N — 3)/2j -times 
continuously differentiable. This result is in agreement with the nonanalytic behaviour 
of the exact solution for the density of states of the mean- field spherical model reported 
in [2]. In other words, the density of states becomes "smoother" with increasing 
number of degrees of freedom, and already for moderate N it will supposedly be 
impossible to observe such finite-system nonanalyticities from noisy experimental or 
numerical data. At first sight one might therefore suspect that the nonanalyticities of 
the entropy are irrelevant for large systems and have no effect in the thermodynamic 
limit, but the following considerations will show that such an assertion is premature. 



4. Thermodynamic limit for the density of states 

In the last years, quite a few articles have been published on the relation between 
nonanalyticities of the entropy or the free energy in the thermodynamic limit and 
topology changes in configuration space (see [2TJ [10] for a review). The object of 
study in this approach is the family {A4v}veTR, with 

My = e Tat I Viq) «C Nv} , (33) 

i.e., the subsets of all points q from configuration space r^v for which the potential 
energy per particle, V{q)/N, is equal to or smaller than a given value v. For several 
models, topology changes of were studied under variation of the parameter v, 
finding that in many cases the occurrence of a phase transition is signalled by a 
signature in the Euler characteristic of My (which is a topological invariant) p2l I23j . 
Then, in a recent letter |7j, Franzosi and Pettini proved that, for a certain class of 
short-range models, a topology change in {A4v}ve'R at v — vt is necessary for a phase 
transition to take place at a transition potential energy vt- This result demonstrates 
that a relation between nonanalyticities of the entropy or the free energy in the 
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thermodynamic limit and topology changes in configuration space does exist at least 
for this class of short-range models. 

It is a central proposition of Morse theory [ini [T^ that, for potentials V having 
the Morse property, each topology change of A^^, at some value v = Vc corresponds to 
one or several critical points Qc of V with critical value Vc = V{qc)/N. More precisely, 
if y is a Morse function and if we know all the critical points of V and their critical 
indices, the handle decomposition theorem [19j asserts that we have all the information 
readily available to specify the topology changes that occur in the family {Mv}v<£-r 
of configuration space subsets. As a consequence, the study of topology changes 
within {Mv}v£Wi. and the study of saddle points of V are closely related approaches. 
From the proven relation between phase transitions and topology changes [7] we can 
therefore conclude that a relation between saddle points of V and nonanalyticities of 
the entropy exists not only for finite systems (as worked out in Sec.[3|), but also in the 
thermodynamic limit of infinite system size. 

In light of our finite-system results summarized in Proposition [T] and Theorem [1] 
this is a remarkable and somewhat surprising observation: Despite their decreasing 
strength for large system sizes N, nonanalyticities of the finite-system entropy appear 
to be related to their infinite-system counterparts in some way. This relation, however, 
is not one-to-one: critical points of V are necessary, but by no means sufficient for 
a phase transition to occur. For several models studied, the number of critical levels 
Vc of the potential V was found to increase unboundedly with the number N of 
degrees of freedom of the system, and the levels become dense on some interval in 
the thermodynamic limit [221 ESI [l3] . Therefore, most of the critical points of V are 
not related to the phase transition. From these observations, the question arises how, 
and under which conditions, nonanalyticities of the finite-system entropy may give 
rise to a phase transition in the thermodynamic limit. This issue will be addressed in 
the present section. 

To this purpose, we investigate the density of states in a small interval 
{vq — e,VQ + e) around some value Vq of the potential energy. In order to quantify 
the contribution to the entropy caused by the critical points, we split the density of 
states into two terms, 

n^riv) = ATi^) + BT'iv)- (34) 

B^''^ contains the nonanalytic contributions, as specified in Theorem [U from all 
critical points in the e- neighbourhood of vq, 

BT(v)^ E (35) 

{v^:\v^-Vo\<e} {q^:V (q^) /N=v^} 

written as a sum over critical values Vc of V inside the neighbourhood and over critical 
points qdvc) corresponding to the respective critical value. Then a smooth function 
A"j^''^ can be chosen such that the e-density of states ft^j^''^ from coincides with the 
exact density of states fijv inside the interval {vq — e,VQ + e). In the following we will 
compute a bound on B^'"^ in order to examine under which conditions the nonanalytic 
terms stemming from the saddle points may yield a non-vanishing contribution to 
n^'*^, and therefore to ^n, in the thermodynamic limit, thereby possibly inducing a 
phase transition. 

In order to perform the thermodynamic limit, we consider the e-entropy per degree 
of freedom, 

s^"^%v) = lim 1 In [A^j^'^v) + B^^iv)] , (36) 
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since, at least for a system with short-range interactions, we may expect this quantity 
to exist [53]. Within the interval {vq — e,Vo + e), s""'^ coincides with the exact entropy 

s{v) = lim — InriAr(u). (37) 

It is important to note that In B^'"^ has in general a non-zero thermodynamic 
limit: For fixed /, we find that the sumniands in Eq. ([50)1 are exponential in N, 



permitting i?^'*^ to give a finite contribution to Now we would like to deduce 
the (non)analyticity of the entropy s from the properties of the thermodynamic limit 
expressions 

a''-^(«)= hm l-\nA"^'\v), (39a) 

N^oo I\ 

b^o^^{v)=lm, ^\nBl"-\v). {39b) 

N—*oo iV 

To this end it is instructive to rewrite Eq. ([36]) as 

s''°''iv)=insiji{a^'>-'{v),b^'"'{v)}. (40) 

The above equation holds unless A^'*^ and B^ '^ are very closely related (such as 
^jv'^ = —B^n''^) which is highly non-generic and thus neglected. As a consequence 
of Eq. (|40|) . we cannot expect to draw any conclusions on s'"°''^ from the knowledge 
of one single function a""'*^ or 6"°''^, but only from the interplay between both. If, for 
example, we find a nonanalyticity in h""'"^ , this nonanalytic behaviour may be visible in 
gfcc jjj one case, but it may simply be overruled by a larger a""'*^ in another instance. 
Furthermore, a nonanalyticity in s'"°''^ may arise from a crossover between a^"'"^ and 
6"°'*^ when the dominant part in the maximum changes from a'"°''^ to 5"°'*^ (or vice 
versa)!]] 

The "density splitting" in [2D], similar to our Eq. raises the hope that 

one may expect the analytic part in A^^'"^ to converge, inside the interval (wq — 
e,fo + e), uniformly to an analytic function for well-behaved short-range potentials. 
Unfortunately we were not able to prove such a result. Furthermore it appears 
plausible from the discussion in the preceding paragraphs that, at least for short- 
range models as covered by the theorem in [7], nonanalyticities of the function 6""''^, 
stemming from the nonanalytic contributions of the critical points of the potential V ^ 
are crucial for the occurrence of a phase transition in the thermodynamic limit. As a 
consequence of this reasoning, we investigate in more detail the thermodynamic limit 
A'^ ^ oo of B]^''^ in order to find out under which conditions it may give a nonanalytic 
contribution to the entropy. 

When performing this limit, we will for simplicity restrict the discussion to 
the subsequence with A^ = 1 (mod 4) and to the leading order |/| = (other 

§ A similar splitting of the density of states into two terms, one stemming from the vicinities of the 
critical points, the other one containing the rest, has been used in |20J. 

II This possibility falsifies the reasoning in I20| saying that, since a first term from volume splitting 
is smooth, a nonanalyticity must be due to the second term. 
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subsequences and higher orders can be treated in a similar way). Now the functional 
form of il'^ can be inserted into the expression yielding 



Here, 



Nr{N/2) \ 



E, N(Ar-2)/2 
[v-VcY 

{vc:va — f.<Vc<va} 



J2 Ji<ic)- E ■^(5^) 



+ E i-c~vr-'^/'\ E Ale)- E 



(41) 



h{v,) = {q, I V{q,)/N = «e A fc(<Zc) = ^ (mod 4)} 



(42) 



(43) 



is the set of all critical points on the critical level which have index k = i (mod 4), 
k{qc) denotes the index of the critical point gc, and 

Nt{v')=Y.A<lc)5{V{qc)/N -v')5tM<i.)rao,iA. £=0,1,2,3, (44) 



are the distribution functions of the critical points with index k — £ (mod 4), 
weighted by their Jacobian determinant J. We expect these Jacobian distributions to 
have a smooth thermodynamic limit in the sense that their integral means 

'^^\v'Afi{v') = - E -^(^c) (45) 



qceQeilv,v+e]) 



have smooth limits relative to the number of critical points in the interval [v,v + e]. 
This implies that 

Je{v) = ^ — (46) 

9ceQf(b,t'+c]) 

has a thermodynamic limit of the form 

J,(w) =exp[7Vj,(z;)] (47) 

for small e > 0. In this case we may express JVg in terms of and the density q of 
critical points with index £ mod 4, 



qaeQi{[v,v+e]) 



yielding 
where 



A^.= E 1' 

qoeQeiR) 



£ = 0,1,2,3, 



(48) 

(49) 
(50) 
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is the total number of critical points with index i mod 4. Generically this number 
grows exponentially 

Ni = exp{Nni), (51) 

for large N, although for specific examples (see the spherical model in Sec. [5]) the 
growth pattern may be different. Collecting all the above information, we obtain an 
upper bound for the contribution of the critical points of V to the density of states, 

BT'iv) < |^^max|7Vi£"^dT;'(«-z;')(^-^)/^ci(«')^i(t''), 

m / dv'{v-v'Y^-^y^cs{v')Mv'), 

Jva~<L 

N2 / Av'iv' -v)^''-^'^'^C2{v')J2{v'), 

Ni dv'iv' -vY^~^^^^c4v')Mv')\. (52) 

Jvo ) 

The fact that the Afi enter Eq. (|41]) as a difference may in very smooth setups lead to 
a B^''^ that is smaller than the right hand side of Eq. ([52|) but we will not need this 
for the arguments presented hereljl Inserting Eq. (|47|) in the above integrals leads to 
Laplace integrals of the form 

In \v - v'\ + je{v') 



f{v')exp{Ne 



2 



(53) 



with some function /(«'). In the thermodynamic limit, these integrals can be evaluated 
by Laplace's method, yielding the maximum of i In \v — v'\ + je{v') within the domain 
of integration. Hence, in the thermodynamic limit (j396p we obtain 



'(v) ^ y/2ne + maxl max ni H — ln(u — v') + je(v') , (54) 

te{0,2},v-e<v' <v 



max 

^6{l,3}, v<v'<v+e 



+ i ln(w' - v) + je{v') 



max [ne+jeiv')] (55) 

££{0,1,2,3}, |t)-u'|<e 

for the contribution of the saddle points in the interval (vq — e,vo + e). 

The implications of this bound on b'"°''^ are the following: From Eq. (HO)) , we know 
that only the larger one of the terms a"°''^ and b'"°''^ contributes to the entropy within 
the interval (wq — e, + e) . Now, if the max-term in (j55p is finite, we can always choose 
e sufficiently small such that the bound ([55|l for b'""''^ is smaller than a^°''^. Therefore, 
the entropy s in the interval {vq — e, + e) is exclusively determined by a""'*^, which we 
assume to be a smooth function. This argument can be repeated for any value of vq , 
thereby demonstrating that the saddle point contributions in the 6-term are harmless 
in the sense that they do not lead to nonanalyticities in the entropy. This reasoning, 
however, relies on the assumption that the max-term in (|55p is finite, and it breaks 
down in case rii or jc diverge for any £ £ {0, 1, 2, 3}. We summarize these observations 
in the following theorem. 

Theorem 2. The saddle point contribution b^°''^{v) cannot induce a phase transition 
at any potential energy in the interval (vq — e, + e) if 

^ The reasoning in on this issue is erroneous. 
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(i) the number of critical points is bounded by exp{CN) for some C > and 

(ii) the Jacobian densities have a thermodynamic limit of the form with 
J, <ooV^e{0,l,2,3} 

inside the given interval. 

Proof. For the leading order, the proof was already sketched above. The leading 
order describes the 'sizes' of the saddle points, whereas higher orders account for their 
detailed shapes. By the geometrical meaning of the term B^'^ it seems implausible 
for a singularity to be caused by the shape and not by the size of the saddle points. 
We do not attempt to give a full mathematical proof for this intuition (abusing the 
term 'theorem'). The result will be corroborated by the examples in Sees. [5] and [7] 
where it is shown that a phase transition is caused by a divergent Jacobian density 
ji{v) at the potential energy of the phase transition v = vt. □ 

Geometrically, a divergent Jacobian density signals that the corresponding saddle 
points become in some sense "asymptotically flat" in the thermodynamic limit. To 
illustrate the implications of these considerations, we calculate in the following sections 
the Jacobian densities je for various spin models in order to check whether the 
requirements of Theorem [2] are met. For the spherical model with nearest-neighbour 
interactions, our results will help us to explain the puzzling observations made by 
Risau-Gusman et al. [13] that, for this model, the topological signatures observed do 
not coincide with the phase transition energy. For the mean-field XY model and 
the mean-field fc-trigonometric model, we will show that a phase transition occurs 
precisely at the value of v for which condition (jn} of Theorem [2] breaks down. 



5. Spherical model with nearest-neighbour coupling 

This model was introduced by Berlin and Kac |25j as an exactly solvable caricature 
of the Ising model of a ferromagnet. Its configuration space Fa? is an A^-sphere with 
radius \/N + 10 and the potential is 

1 ^ 

Vsph : r^v -> IR, ^ ^2 ^ K,jqiqj. (56) 

The degrees of freedom qi of the model are placed on a d-dimensional cubic lattice 
with -I- 1 = L'^ lattice sites and periodic boundary conditions. We consider a 
coupling matrix K with elements Kij = 1 when i and j are nearest-neighbouring sites 
on the lattice, and Kij = otherwise. The spherical model with nearest-neighbour 
coupling is solvable in the thermodynamic limit for arbitrary d, and a phase transition 
from a ferromagnetic phase at low temperatures to a paramagnetic phase at high 
temperatures (or energies) occurs for all c? ^ 3 [53]. Since the coupling matrix K is 
symmetric, it can be diagonalized by a constant orthogonal transformation, and we 
can rewrite the potential in the form 

1 ^ 

Vsph : Fat ^ R, ^ ~ 2 ^ '^*^*^' ^^'^^ 

i=Q 

+ We shift N to an unconventional A'' + 1 to match notation with the previous sections. 
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where (i = 0, . . . , TV) are the eigenvalues ofK. As discussed in [13], these eigenvalues 
can be written as 



X{pi, . . . ,pd) = 2^ cos f^Y^ j , 



p, = 0,...,L-l, (58) 



e=i 

and we consider the Xi in ([57|) as the values oi X{pi, . . . ,pd), ordered from the largest 
to the smallest. Several of these eigenvalues are degenerate, and we define the spectral 
density of K (as a function of the potential energy per particle v) as 

1 

<^jv(^) jv + 1 ^ S{X{pi,...,pd)-v). (59) 

pi,...,Pd=0 

In two recent publications, Risau-Gusman et al. [131 114] reported a study of the 
topology of the subsets Aiy C r^r [as defined in ((33|) ] of the spherical model with 
nearest-neighbour coupling. Computing the deformation retract of the M^, they 
found that their topology can be characterized by means of the spectral density (j59p . 
Therefore a signature of the topology changes in the limit of large system size should 
be visible in the thermodynamic limit TV + 1 = L'^ ^ oo of the integral mean of the 
spectral density ctjv- It was shown in [14j that this mean converges to 

1 f°° 

lim CTjv(w) - / dxcos(2a;v)[Jo(a;)]'', (60) 

where Jq denotes the Bessel function of order zero. Depending on the spatial dimension 
d, this limit function has one or several nonanalytic points. However, none of these 
coincides with the phase transition energy (see Figure 1 of Reference for a 
plot of these functions for various values of d). Therefore no signature of the phase 
transition of the spherical model is seen in the topological quantities studied. 

Inspired by the analysis in [HI [14], we study the subsets Al^ as defined in ([33)l 
by means of Morse theory. Although the potential Vsph of the spherical model in ([57| 
is not a Morse function, it can be transformed into one by adding — in the spirit of 
Eq. dU) — perturbations hiX^ with small hi € R, 

1 ^ 

V^phtTAT^R, --Y^X^xl (61) 

i=0 

where Ai = Ai + Zi^ (i = 0, . . . , N). For convenience, the hi are chosen such that Xi > Xj 
for i < j. For this modified potential we want to calculate the Jacobian densities J^i 
as given in (|46p when performing simultaneously the zero-field limit h and the 
thermodynamic limit N — > oo. To this end, we consider hi as iV-dependent quantities, 
vanishing in the thermodynamic limit in a suitable way. 

As a first step, we determine the critical points and critical indices of Kiph, using 
the method of Lagrange multipliers. In the Lagrange function 

N N 

£ ^ ^ Kx', +^,[j2x^-{N + 1)] (62) 

1 = i=0 

the second term on the right hand side takes into account the spherical shape of the 
configuration space of the model, and fi is the corresponding Lagrange multiplier. The 
critical points of Vsph are obtained as the solutions of 

d£ 

— = -X,Xi + 2fix, = \/i = 0,...,N. (63) 

OXi 
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Simple algebra shows that 2(iV + 1) critical points exist which can be written as 

xi^''^ ^±VWTTek, k = 0,...,N, (64) 

where is the unit vector in fc-direction. Inserting this expression into the potential 
(|6ip . we obtain the corresponding critical values 

y,^YE2hi^^^^(^l + N~'), k^O,...,N, (65) 

and we note that precisely two critical points are located at each critical value. The 
density of critical points is thus linked to the spectral density crjv, and in the limit 
N + 1 = L"^ ^ oo the integral mean of the density of critical points is 

c{v) = hm (Ttv f~^) ~ ~ I dxcos(xu)[Jo(a;)]'', (66) 



similar to Eq. ((60|) . 

Now we fix a certain value of fc G {0, . . . , N}. A good set of coordinates in the 
vicinity of the critical values Xc is given by the Xi with i ^ k. Solving the spherical 
constraint 

N 

Y,xj-{N+l) = (67) 

for Xk and inserting the resulting expression into the potential (j6ip . we obtain 

V.^x) ^-\Y.^K- \k)xl - ^^Afe. (68) 

The potential is quadratic in the Xi and leads to a diagonal Hessian at the critical 
points with matrix elements 

^^,{x^t'''^) = (Afc - \{)5,, for i ^ k. (69) 

Due to the ordering of the , we can read off from this expression that the index of 
a critical point Xc equals k. Now we can write the Jacobian determinant in terms 



of the critical values Vi ([65]) , 

J(4^^)) = |det[^(xP))/2]r^/^=n ^^'''-expjlj: In 1.(70) 



Expressing the Jacobian determinant as a function of the potential (per particle) v 
allows us to formulate the thermodynamic limit in terms of the density c{v) of critical 
points (|66l) . 

^ r+d 

J2 dvciv), (71) 

where d is the spatial dimension of the lattice on which the spherical model is defined. 
The range of integration [— d, +d] is determined by the minimum and the maximum 
values of the potential energy per particle. After some manipulations of the integral 
obtained from ((70)) with the substitution ((7T1) . we arrive at the Jacobian determinant 
in the thermodynamic limit as a function of u, 

J(v) = e^p[-— dx[Mx)ff{v,x)j. (72) 
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The function / in this expression is given by 

xf{v, x) = sm{dx) In (d^ - u^) + sin(a;i;) [Ci{x{d + v)) - Ci{x{d - v))] 

- co8{xv) [Si{x{d + v)) + Si{x{d - v))] , (73) 

where Ci and Si are the cosine and sine integral functions. f{v,x) is an analytic 
function for all accessible values v £ {~d,d). As a consequence, also S{v) is analytic 
and, in particular, continuous for all v. 

Observing that for the spherical model both, the index of a critical point Xc 
as well as the Jacobian determinant J(xc^''^), are uniquely determined by the critical 
values Vk = Vsphixc^'^"' ) / N which decrease monotonously with increasing j, we find 

Ji^J and Q = c, ^ = 0,1, 2, 3. (74) 

We thus may apply Theorem[2] (with any C > because Ne — > (A^+l)/2) and conclude 
that for the spherical model with nearest-neighbour coupling the contributions from 
the critical points of the potential (or, equivalently, from the topology changes of My ) 
to the entropy are negligible in the thermodynamic limit. 

At first sight this finding seems to be in conflict with the previous discussion where 
the relevance of the 5-term in ([40]) for the occurrence of a phase transition in systems 
with short-range interactions was emphasized. This apparent contradiction can be 
understood by noting that the spherical model is not genuinely short-range. Although 
the interaction term in (|56p is restricted to pairs of nearest neighbours, the spherical 
configuration space renders the interaction effectively long-range. For long-range 
systems, however, as discussed in [27], the connection between phase transitions and 
configuration space topology is not valid in general and we cannot assume the function 
a""'*^ in (j40p to give a smooth contribution to the entropy in the thermodynamic limit. 
Similar to the observations made in |27| for a different long-range model, we suspect 
the entropy s{v, m) as a function of the potential energy v and the magnetization m to 
be a nonconcave function, such that the nonanalyticity in s{v) = max^ s(v,m) arises 
from the maximization over m. An analytic calculation of the entropy s{v,m) of the 
spherical model with nearest-neighbour coupling has been reported in '28], but this 
analysis omits precisely the region in the (u, m)-plane for which a nonconcavity of the 
entropy might be expected. 

Note that, in the presence of long-range interactions, we can not a priori 
guarantee a phase transition to be triggered by saddle points of V (or topology 
changes of Aiv), since the conditions of the Franzosi-Pettini- Theorem are not met. 
Nonetheless, topology changes may be at the origin of a phase transition even in 
long-range systems, and we present two such examples in the following sections. 



6. Mean-field XY model 



The mean-field XY model is a system of N plane rotators described by angular 
variables qi G [0, 27r) {i = 1,...,N). Each rotator is coupled to each other with 
equal strength K > 0, where the interactions are described by the potential 

K ^ ^ 
VxY ■■[0,2Trf ^R, ^[l_cos(g,-g,)]-/i^cos<7,. (75) 

— l i—1 

The rotators are subject to an external magnetic field of strength /i G IR which 
energetically favours orientations qi ~ n for h < 0, respectively w (mod 27r) 
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for h > 0. In the limit N ^ oo and — > 0, the system has a continuous phase 
transition. 

The critical points of this model and their indices have been analyzed by Casetti 
et al. [221 for an arbitrary number iV G N of rotators. Critical points of VxYio) 



were found to occur for q 



e {0, tt} , i. e., for {qi, . . . , q^) with all components 



qi {i — l,...,iV) being either or tt. Hence, in accordance with condition ^ of 
Theorem [21 the number of critical points is growing exponentially with N . The 
elements — d^V/{dqidqj) of the Hessian are also reported in [22] and, evaluated 
at a critical point, the diagonal elements of the Hessian can be written as 



K 



N 



N ^ 

fc=i 



cos qk + n] cos 9i - = 



K 1- 



cosqt-^, (76) 



where n-^ is the number of tt's in the sequence (gi, . . . , gjv) € {0, tt}^. It has been 
shown in Appendix A.l of [22] that, in the thermodynamic limit N ^ oo, the 
contribution of the off-diagonal elements to the eigenvalues of Sj vanishes. Therefore 



we can write the Jacobian determinant 



J(g(^)) = det[fi(q('=))/2] 



I at a critical point q^'^^ as 

^1/2 ^ 



4=1 



^11- 



-N/2 



valid for large N. Inserting, again from Ref. [22 , an expression of n, 
the potential energy per degree of freedom v, 

n^{v) Ar>i 1 



K 



h ± + - 2Kt 



N 2K 

we obtain the Jacobian determinant as a function of v, 



-1/2 

(77) 

as a function of 
(78) 



J{v) 



N/4 



+ h? - 2Kv 
valid for large N. In the zero-field limit h 

4 



this expression reduces to 



J(.v) 



N/4 



(79) 



(80) 



K{K-2v) 

Similar to the calculation for the spherical model in Sec. [5l the indices of the critical 
points increase monotonously with their respective critical values. Hence, we can 
conclude that Ji — J (t = 1, 2, 3) and, making use of Eq. (|47p . we obtain 



ln2- ^\xi\K(K- 



2^^)], 



£ = 0,1,2,3. 



(81) 



2 4 

The graph of this function is plotted in Fig. [1] As is easily recognized, ji(v) is 
singular only at u = A'/2, which is precisely the value of the potential energy per 
degree of freedom at which a phase transition occurs for the mean-field XY model 
with zero external field in the thermodynamic limit. This result nicely illustrates the 
content of Theorem [2J no phase transition is triggered by the critical points of the 
potential with a finite Jacobian density. Only those critical points for which becomes 
singular may induce a phase transition, and this is what apparently happens for the 
mean-field XY model with /i = 0. Note that, for h and in the thermodynamic 
limit, entropic reasons prevent the system from visiting states with potential energy 
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V = [K"^ + h?') / {2K) for which the Jacobian determinant is singular. As a consequence, 
in accordance with the known thermodynamic behaviour of the mean-field XY model, 
no phase transition occurs in the presence of a non-zero external field h. 

Note that, strictly speaking, the analysis of isolated critical points as presented in 
Sees. |3] and [4] cannot be applied to the mean-field XY model at the maximum critical 
level V = [K^ + h?)l{2K): Precisely at this level the potential Vxy is not a proper 
Morse function, as its critical points are degenerate for arbitrary TV and form a critical 
manifold. However, this should not lead to serious problems for the interpretation of 
our result: Excluding the trouble-making maximum critical level from our analysis, we 
can monitor the behaviour of the Jacobian determinant for the second largest critical 
level. For any finite N , the critical points on this level are isolated and hence our 
analysis is applicable. In the thermodynamic limit, this second largest critical level 
approaches the problematic maximum one arbitrarily close and the corresponding 
Jacobian determinant diverges. Note, however, that a divergent Jacobian density 
need not necessarily be related to the break-down of the Morse property: One might 
as well imagine a system where all critical points are isolated for all finite TV, but for 
which nonetheless the Jacobian density diverges at some point in the thermodynamic 
limit TV — > oo. 

7. Mean-field fc-trigonometric model 

Another spin model, similar to the one discussed in the previous section, is the mean- 
field fc-trigonometric model. A simplified version of this model was introduced in 
as a model of a simple liquid. Like the XY model, the system consists of TV plane 
rotators described by angular variables qi € [0, 2tt) [i = 1, . . . , TV). The Hamiltonian 
function of the mean-field fc-trigonometric model is 

N 

Vk : [0, 2^)^ R, ATV^"^ ^ [1 - cos(g,, + . . . + J] , (82) 

where q= {qi, . . . , qN)- The constant A determines the coupling strength between the 
rotators, each interacting with each other at equal strength. Both, thermodynamic 




Figure 1. Plot of the graph of the Jacobian density jf as a function of the 
potential energy v for the mean- field XY model with coupling constant K = 1. 
In agreement with Theorem [J] the transition potential energy vt = 1/2 of this 
model coincides with the singularity of ji . 
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behaviour and critical points of the potential have been studied by Angelani et al. 
[3D], and the following results are taken from this reference. 

The parameter fc g N in the Hamiltonian function crucially determines the 
thermodynamic behaviour of the model: In the limit TV — > oo, the system shows 
a discontinuous phase transition for fc ^ 3, a continuous one for fc = 2, and no phase 
transition for fc = 1. Results on the critical points of Vk and on the Hessian at these 
points can be assembled to obtain the Jacobian density ji as a function of the potential 
energy per particle v, 



je{v) = iln 



2 

Afc 



1 



V 

A 



l/k-l 



£ = 0,1,2,3. 



(83) 



The graph of this function is plotted in Fig. [2] for k = 1,2,3,4. The calculation is 
analogous to the one for the mean-field XY model in the previous section and we 
hence skip the details. 
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Figure 2. Plot of the graph of the Jacobian density ji as a function of the 
potential energy v (in units of the coupling constant A) for the mean-field k- 
trigonometric model with k = 1,2,3,4. In agreement with Theorem |2] the 
transition potential energy vt = A o{ this model coincides with the singularity of 
jl for all values k ^ 2 for which a phase transition takes place. 

This is another result nicely illustrating the content of Theorem [2j The function 
jt is constant (and therefore bounded above) for fc = 1 , in agreement with the fact that 
no phase transition occurs in this case. For fc ^ 2, however, the Jacobian density jt{v) 
shows a divergence at w = A, which is precisely the value of the potential energy per 
degree of freedom at which the mean-field fc-trigonometric model undergoes a phase 
transition in the thermodynamic limit. 

8. Summary 

We have analyzed the relation between saddle points of the potential energy V of 
classical A'^-particle systems and the analyticity properties of their thermodynamic 
functions. For finite systems, each saddle point Qc was found to cause a nonanalyticity 
in the entropy sn{v) at the value v = V{qc)/N of the potential energy, and the 
functional form of the nonanalytic term is given in Theorem [T] Since the number 
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of saddle points is expected to grow exponentially with N for generic potentials, 
we arrive at the remarkable conclusion that the finite-system entropy is in general 
a highly non-smooth function. For large N, the order of the nonanalytic term 
increases unboundedly, leading to an increasing differentiability of sn . Considering the 
contribution of very many saddle points becoming dense in the thermodynamic limit, 
we discussed how, despite the "increasing smoothness" of sat, a continuous distribution 
of saddle points with singular Jacobian density may lead to a nonanalyticity in the 
infinite-system entropy. Interpreting our findings in the spirit of the topological 
approach to phase transition, our results show under which conditions topology 
changes of the configuration space subsets as defined in p3|) can be at the origin 
of a phase transition in the thermodynamic limit. 

An application of these findings to the spherical model allowed us to understand 
the puzzling observations of Risau-Gusman et al. [13l [T^ that, for this model, the 
topological signatures observed do not coincide with the phase transition energy. For 
the mean-field XY model and the mean-field fc-trigonometric model we showed that 
a phase transition occurs precisely at the value of the potential energy per degree of 
freedom v for which condition ^ of Theorem [2] breaks down. For these two mean-field 
models, the divergences in the Jacobian densities ji in the thermodynamic limit are 
foreshadowed by the break-down of the Morse property of the potentials already for 
finite systems. Note, however, that this is not necessarily the case and one can as well 
imagine a divergent Jacobian density to emerge in the thermodynamic limit from a 
potential which is a good Morse function for all finite system sizes N. 

The criterion for the absence of a phase transition we have presented in Theorcm[2] 
is to some extend related to what has become known as the topological hypothesis in the 
literature, where a relation between the occurrence of a phase transition and certain 
"pronounced" topology changes of the subsets Mv [Eq. ([33)1 ] in configuration space was 
conjectured (see [21] for a discussion of the different formulations of this hypothesis). 
In fact, for the mean-field XY model and the mean-field fc-trigonometric model in 
the thermodynamic limit, signatures were found in the Euler characteristic of My, 
occurring precisely at the transition potential energy v = Vt, and a relation between 
these signatures and the occurrence of phase transitions was conjectured [121 [53]. 
Our Theorem [21 similar to this conjectured relation, has a topological ingredient as 
well, taking the saddle points of the potential (each of which is related to a topology 
change of the A4v) as a starting point. In addition to this topological ingredient, we 
have demonstrated that the local curvatures at the saddle points are important to 
quantify the effect of these saddle points on thermodynamic quantities. In this way, 
we derive a criterion, in part of topological and in part of geometrical nature, which 
can exclude the occurrence of a phase transition for most of the accessible values of 
the potential energy v, leaving only a few distinct values of v as possible candidates 
for the occurrence of a phase transition. 

Conceptually, our study explores the connection between phase transition theory 
and the study of energy landscapes, a rapidly developing field with applications, among 
others, to clusters, biomolecules, and glass-formers [9]. One might hope to profit 
from the considerable knowledge on the relation of energy landscapes and dynamical 
properties in future work. 
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In this Appendix we will examine the effect of an x-dependent Jacobian J{x) on the 
singular behaviour of the density of states Qat. 

Upon Taylor-expanding the Jacobian, we obtain terms proportional to = 



. These terms result in corrections to Eq. (|T3)) of the form 

Ai^aiJ d^x x'5 {-X^ + Y^-V)Q (r^ - - V') , (A.l) 



where Oq = | det[i3(gc)/2]| Due to the symmetry upon substituting 



we see that only even indices ij contribute. We use 9(r^ — X"^ — 
X^ — Y"^) and express A/ in terms of a generating function 

Aj = ai dyrWdlGip,y) , 
Jo p=o 

G{p,y) = / d^x5{-X'^ +Y'^ -V)6{y- X"^ -Y'^)eyi^{iY^ 



XjPj 



Xj, 

/J" dy5{y- 

(A.2) 
(A.3) 



where we have used the notations = d]^ dl^ ■ ■ ■ d^j^ , dj = d/dpj, and |/| = 
ii + *2 + ••■ + ijv- Now we represent the (5-distribution by their Fourier components. 
For convergence we introduce a small e > and obtain 

1 



G{p,y) 



(2^)2 



ds I dt I d^cc 



N 



X exp 



s{-X^ + Y^ -V) + (it + e){y - X^ - Y^) + iJ2 ^jPj - + ^' 
The argument of the exponential is quadratic in x and can be rewritten as 



(A.4) 



k 

E 

i=i 

N 

E 

j=k+i 



-{e + is + it) I X 



ipj 



2(e + is + it) J 4(e + is + it) 



^(e — is + it) ( Xj — 



\pj 



2(e - is + it) / 4(e - is + it) 



-isV +{it + e)y-e{s'^ +t^). (A.5) 
From this expression we read off the value of the Gaussian integral and obtain 



Gip,y) 



(2^)2 



dt 



exp 



-N-k 



„ „ (A.6) 

ye + IS + it y'e — is + it 

using the same definition of the complex square root as in the main text. The 
generating function factorizes in the variables Pj. We use the shorthands Ii = X]^=i *j 
and I2 = X^jL/c+i *j have Ii + I2 = \I\) and obtain 



ai 



(27r)24l^l/2n;Li(*,/2)! 

exp [-isy + ity - e{s'^ + t^)] 



dye"^ ds dt 



I : rrfc+Zl , : tt 

ye + IS + it ye — is + it 



—N-k+I-i ■ 



(A.7) 
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Now we observe that the above double integral in s and t factorizes after changing 
variables to u ~ {t + s)/2 and w = (i — s)/2. We define the function 

1 r , 1 



F,{z) = ^ du ^ e--^^- , (A.8) 

■^''^ J~oc ve + iu 

drop the factor e'^^ which has no effect on the finite integral over y in the limit e ^ 0, 
and replace e by 2e, yielding 



2a,(V2)^/^n,ti», 
8l^l/^nf=i(*./2)! 



dy Fk+i, (y - V)FN-k+i, (y + F) . ( A.9) 



We have to evaluate Fq{z) for real values of z and non- negative integers j. Obviously, 
F^{z)^5{z). 

For positive, even j = 2n the integrand has a pole at u = ie of order n. If 2: < 
we may evaluate the integral by closing in the lower half plane. We miss the pole and 
F2n vanishes. If z > we may evaluate the integral by the residue theorem closing in 
the upper half plane. Expanding the exponential gives — 0(z)z"~^/(n — 1)!. 

In the case j — 1 the integrand has a cut from u = ie to ioo. Again we obtain a 
non- vanishing result only in the case z > where we may deform the contour to follow 
the cut. We let e approach zero and obtain Fi{z) = -R-^Qiz) due'"Viv/H- This 
leads to a Gamma- function evaluated at 1/2 and we have Fi = {■kz)-'^/^Q{z). 

The case of odd j = 2n-\-l can be reduced to the case j = 1 by integration by parts. 
We may drop the factor e~^'^" in this case and obtain F-2n+i{z) = F2„_i(z) 2z/ (2n— 1). 
We iterate this equation and use the above result for j = 1 to obtain a result analogous 
to the case of even j . This allows us to write the result in compact form 

Fo(z) = J(z), F,{z) ^ ^^^<d{z) ifj^l. (A.IO) 

Now we have to distinguish the cases of a maximum, a minimum, or a proper saddle 
point. In the first case we have k — Ii = Q which reduces Fk+i^ [y — V) to 5{y — V). 
We only obtain a non- vanishing result if < < r^, and collecting all factors gives 

A, — ^-^^^ l-(^+l^l-2)/2e(F)e(r2 - V) (A.U) 

2mr[(A^ + |/|)/2]nf=i(*,/2)! 

for A: = 0. In the case of a minimum we have k = N, I2 = with a non-vanishing 
result only for — < V < 0. Upon V —V we reproduce — as it should be — the 
formula of a maximum 

A, — — ^^'=' '- (-y)(~+i^i-2)/2e(-y)e(r2 + v) (A.12) 



2l^lr[(A^+|/|)/2]n;=ife/2)! 

for k = N . In both cases we confirm the statements of Proposition [T] and Theorem [T] 
The case of a proper saddle yields 



2aj{n/2r'^Y{Uh 



m/^V [{N + h)l2] T[{N-k + h)l2] nf=i(^,/2) 



(A.13) 



dy (y - T/)('=+^i-2)/2(y + y^(N-k+I,-2)/2 (^^^4) 
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for < k < N . This expression entails the integrals I± from Eq. (22) in a symmetrized 
form. In fact, a transformation y = V{2z — 1) allows us to express A/ in terms of I± 
with k replaced by fc + /i and N replaced by iV + |/| (denoted by an upper index), 



.1 



2WT[iN + h)/2]T[{N-k + h)/2]UU{^,/2)\ ^ 

Now we can follow the discussion in the main text substituting {k, N) by (k+Ii, N+\I\) 
when necessary. We observe that, with Ii and |/| being even integers, higher order 
corrections do not mix the cases in the derivation of Proposition [1] Collecting all 
factors gives the desired result. 
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